Freezing transition in a direction-free kinetically constrained 
model on a finite dimensional lattice 
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Abstract. - We propose a direction-free kinetically constrained model on the square lattice 
and show that this model exhibits a freezing transition at a non-trivial density. By numerical 
experiments, we also conjecture that critical exponents characterizing the singular behaviour near 
the transition point are different from those of previously known freezing transitions in finite 
dimensions. 
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Introduction. — Athermal particles with repulsive 
interactions exhibit rigidity with an amorphous structure 
when the density of the particles is higher than a criti- 
cal value. Such a phenomenon is referred to as a jam- 
ming transition, and elucidation of the nature of jamming 
transitions is a challenging problem in statistical physics 
[IHS]. Thus far, the replica theory in equilibrium sta- 
tistical mechanics has provided an insight into jamming 
transitions as well as glass transitions [Hllin]- As a dif- 
ferent approach, kinetically constrained models (KCMs), 
. which were investigated with a physical picture that glassy 
■ dynamics is purely kinematic [llj , have been considered 
for understanding jamming transitions |12j . A notable 
' property of KCMs is the absence of singularities in the 
. equilibrium properties. Although KCMs are simple math- 
' ematical models that lack several properties of realistic 
. materials, recent experimental studies have attempted to 
' reveal the relationship between KCMs and glass-forming 
materials P^fTi) . 

Let us review theoretical studies on KCMs in brief. The 
first rigorous proof for the existence of a freezing transition 
was presented for the Fredrickson- Andersen model and the 
Kob- Andersen model on a Bethe lattice JlSJ . In this case, 
a freezing transition means a transition to a freezing phase 
where an infinite number of particles are at rest as a re- 
sult of blocking by other particles. However, for these 
models in finite dimensions, it has been shown that there 
are no true phase transitions [16l[T7], even though the re- 
sults of numerical simulations suggested the existence of 



transition. Then, a remarkable discovery was made that 
certain KCMs exhibit a freezing transition in two dimen- 
sions [T2 jlT8H20] . Here, it should be noted that the con- 
straint conditions of these models depend on directions. 
This characteristic property of KCMs distinguishes them 
from the more familiar direction- free models such as the 
Ising models. 

The basic motivation in the present study is to estab- 
lish the idea of universality classes for freezing transitions 
observed in finite-dimensional KCMs. However, this prob- 
lem cannot be considered directly, because such KCMs are 
scarcely found. Hence, we first need to collect examples 
of KCMs that exhibit a freezing transition in finite dimen- 
sions. In particular, because all the examples known thus 
far belong to the same universality class, which we refer to 
as the Toninelli-Biroli- Fisher (TBF) class, it is a signifi- 
cant challenge to present a KCM that does not belong to 
the TBF class, and exhibits a freezing transition despite 
this. We attempt to solve this problem. 

In particular, we focus on direction-free KCMs, which 
may be conjectured as qualitatively different from KCMs 
in the TBF class because of the symmetry property re- 
lated to directions. Direction-free models also present an 
advantage in that their mean-field analysis may be easier. 
Mean-field analysis acts as a starting point for a theory 
on determining the universality, as we have learned from 
the renormalization group analysis for equilibrium criti- 
cal phenomena. Therefore, the discovery of a direction- 
free KCM that exhibits a freezing transition might be 
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Fig. 1: Two configurations satisfying XljgB "^i ~ 2. In the 
configuration to the left, the particle at site i is constrained, 
because fj and fji axe at most 1. In the configuration to the 
right, the particle at site i is not constrained when there are 
no particles at the sites marked with the star symbol. The 
particle at site i becomes constrained when a particle is placed 
at the sites marked with the star symbol. 



useful for further studies on the universality. With this 
background, we propose a direction-free two-dimensional 
KCM for which the existence of a freezing transition can 
be proved rigorously. We also conjecture by numerical ex- 
periments that the proposed model does not belong to the 
TBF class. 

Model. — Let A be a square lattice consisting of sites 
i = {ij:,iy), where ix and iy are integers satisfying 1 < ix < 
L and 1 < iy < M. We assume M < L without loss of 
generality. We define the occupation variable cr^ G {0, 1} 
at each site i G A, and we assume a simple Hamiltonian 



(1) 



where we express cr = (cri)igA collectively. Concretely, 
(Ti = 1 represents a state that a particle occupies site i, 
while there is no particle at site i when (7^ = 0. In the 
argument below, we assume that particles are filled outside 
the system unless otherwise mentioned explicitly. 

We consider a Markov stochastic process with a tran- 
sition ratio R{cr cr') for cr cr' . Let P{cr,t) be the 
probability distribution of a particle configuration cr at 
time t. The time evolution of P{cr,t) obeys 

dtP{(T,t) = [R{ct' cr)P{(T',t) - R{cr cr')P{(T,t)]. 

In this work, we study a constrained Glauber dynamics, 
where the transition ratio R{cr — > cr') is written as 



The term S{cr' , Ficr) represents a rule that the transition is 
caused by a particle creation or annihilation at each site. 
The term r{cr ^ cr') is given as 



r{cr — > cr') = min 

where it satisfies 
r(cr a-') 



1, exp 



H{a)-H{a') 



r{a' -)■ (t) 



exp 



Hia')-Hia) 
T 



(5) 



(6) 



Finally, with regard to the function Ci{cr), Ci{cr) — spec- 
ifies a set of configurations for which particle creation and 
annihilation at site i is prohibited, and Ci{cr) — 1, for the 
other configurations. We also assume that Ci{cr) is inde- 
pendent of ffi. 

It should be noted that the canonical distribution with 
the Hamiltonian ([1]) is the stationary solution of Q with 
(jS]) because the transition ratio ([3]) satisfies the detailed 
balance condition. Therefore, there are no equilibrium 
phase transitions in the system. Nevertheless, it has been 
shown that there is another steady state (a freezing state) 
for an appropriately selected Ci(cr), which is different from 
the equilibrium steady state. All examples of such Ci{cr) 
known thus far depend on directions [16l[T8l[20] . In con- 
trast to those in the previous cases, we propose a direction- 
free function Ci{cr) that leads to a freezing state, where, in 
precise terms, a direction- free Ci{cr) should be expressed 
only in terms of B j , which indicates the set of sites con- 
nected to site j in the lattice A. 

Now, we explain our selection of Ci{cr). We first define 



Y,5{a,,Q) 



(7) 



where 5{m, n) represents Kronecker's delta function. For 
site i, fi represents the number of empty sites j in 
such that there are no particles next to any site j. From 
this definition, we find fi = for site i with a particle. 



Then, we set Ci(cr) = 
(ii) E,eB, = 2 and 



0, only when (i) X^jeB; > 3 or 



(8) 



I -^-^ I I satisfying the condition X^jee '^i 

R(cT ^ cr') 2^ (5(cr', F^cT)r{cT cr')C^{cT). (3) ^^^^ gg^j.e, and /. 



We explain the right-hand side of the equation in order. 
First, Fi is the creation and annihilation operator de- 
scribed by 



{Ficr)j = (1 - ai)6ij + aj{l - Si^ 



(4) 



(2) 

Condition (i) is a standard constraint for KCMs, while the 
additional rule (jS]) in condition (ii) may appear unfamiliar. 
We explain rule (ii) through the following examples. 
In figure [U for the particle at site i, there are two cases 

= 2. In the configura- 
for the empty sites j 
and j', are at most 1. As X^jeB; ^(/jj 2) = in this case, 
Ci{cr) = 0. Throughout this work, we say that the particle 
at site i is constrained when Ci{cr) = and Gi = \. On 
the other hand, in the configuration to the right in figure 
[l] when there are no particles at the sites marked with 
the star symbol, fj and fji for the empty sites j and j' 



p-2 



Freezing transition in a direction-free KCM on a finite dimensional lattice 



are 2, and consequently, J2jeBi '^(/ji^) = 2. In this case, 
the particle at site i is not constrained. When at least one 
particle is placed at the sites marked with the star symbol, 
J2jeB ^(/ji 2) = 1 and the particle at site i is constrained 
as the result. Using similar considerations, we also find 
that the function Ci{(T) is independent of (7^. 

In this work, we focus on the relaxation behaviours of 
the system from an initial state, where sites are randomly 
occupied with a probability p. Moreover, as the simplest 
case, we consider the zero temperature limit T — ^ 0. The 
equilibrium state in this case involves no particles. Now, 
let us suppose that a particle configuration at t = con- 
tains a set of particles constrained by only constrained 
particles. Such particles are referred to as frozen parti- 
cles. Then, frozen particles cannot annihilate, and thus 
they can never reach the equilibrium state. We refer to 
this phenomenon as freezing. On the other hand, when no 
frozen particles exist at i = 0, all the particles annihilate, 
and hence, the final state becomes the equilibrium state. 
In this manner, whether a freezing transition occurs is de- 
termined by investigating the existence of frozen particles 
in initial configurations. 

Existence of a freezing transition. In this sec- 
tion, we prove that a freezing transition occurs at some 
value < pc < 1 in the proposed model. In concrete 
terms, we show that there is a density pi > 0, below which 
the probability of finding frozen particles in initial config- 
urations is zero in the thermodynamic limit, while there 
is a density < 1, above which frozen particles can be 
found in initial configurations with probability one in the 
thermodynamic limit. By estimating pi and pu, we may 
conclude that < pi < pc < Pu < 1, which means that 
pi and yOu are a lower bound and an upper bound of the 
transition point, respectively. 

No frozen particles for p < p\. Let (i^'^^)^i be a 
sequence of A^-sites such that (i) a^m = 1, (ii) ji^*^^ — 
I < 3, and (iii) \i^'''> - i^^"^ | > 3 for fc' 7^ fc ± 1, where 
|i — j| = \ix — jx \ + \iy — jy\- For a given configuration cr, 
we define a set of all such sequences, which is denoted by 
'Dn{o')- Here, if frozen particles are present in the initial 
configuration cr, we can always find a sequence of sites in 
T^M/si^)- (Recall that M < L.) Thus, the probability 
Q of finding frozen particles in initial conditions is less 
than Prob(X'j\//3 ^ 0). Here, the latter probability is 

less than (24p)T^, where 24 is the number of sites whose 
distance from one site is at most three. These estimations 
give Q < (24p)-. When p is less than 1/24, Q ^ in 
the thermodynamic limit. We thus find a lower bound as 
pi = 1/24. 

Frozen particles for p > p^. Let i3 be a set of ori- 
ented bonds written as {(4fc, 2/) (4fc -I- 2,2/ ± 1)} or 
{(4fc -2,21 + 1) (4fc, 2l + l± 1)}, where k,l e Z. We 
say that a bond {(n, m) (n + 2, m ± 1)} £ B is occupied 
when three sites (n, m), (n+l,m), and (n-|-l,m±l) are all 
occupied. Now, if there exists an infinite connected clus- 



9000 




0.65 0.66 0.67 0.68 0.69 0.7 0.71 

P 



Fig. 2: Relaxation time tq with L = 10000. 

ter (in the sense of a directed percolation) of the occupied 
bonds in B, we can find frozen particles in the thermo- 
dynamic limit. The problem is similar to the standard 
directed percolation. However, the occupation of bond 
{(4A:,2/) (4fc + 2,2/ + 1)} is not independent of the 
occupation of bond {(4fc, 21) (4fc -I- 2, 2/ — 1)} because 
the occupation of the two sites, (4A:, 2/) and {4k + 1,21), 
influences the two bonds. 

In order to estimate the condition for the existence of 
an infinite connected cluster, we introduce two auxiliary 
variables, and , which take the value with prob- 
ability I — p, and 1, with probability p. By using these 
variables, we may determine an initial configuration cr by 
using the relation at — -I- /3~ — fif f^^ ■ The probability 
p is then related to p as p = 1 — (1 — pf' . Here, instead of 
the bond defined by the occupation variable Ui, we con- 
sider the bond defined by the auxiliary variables [3^, as 
follows. We say that a bond {{n, ni) — > (n + 2,m + 1)} is 
/3-occupied when = = = 1. Like- 

wise, {{n, m) —i' {n + 2,m — 1)} is called /3-occupied when 
/3n,m = /^n+i,m = /^n+i,m-i = 1- ^s /3-occupation for cach 
bond independently occurs with the probability p — p'^, 
there is a critical value pc = p^ — 0.644 • ■ •, above which 
the directed bond percolation occurs [21]. If a bond is 
/3-occupied, the bond is occupied. Thus, an infinite con- 
nected cluster of /3-occupied bonds indicates an infinite 
connected cluster of bonds in the original problem. By set- 
ting pu = 1 — (1 — pc)^, we find that there exists an infinite 
connected cluster (in the sense of a directed percolation) 
of occupied bonds in B iov p > p^. This pu ~ 0.981 • • • is 
an upper bound of p^. 

Transition point and universality class. In this 
section, with the aid of numerical experiments, we ob- 
tain greater lower bounds and smaller upper bounds of 
the transition point than the rigorous estimations given in 
the previous section. We also present evidence that the 
singular behaviour near the transition point is different 
from that for the TBF class. 

Greater lower hounds. Lower bounds of the transition 
point are numerically obtained by checking whether initial 
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Fig. 3: Pl(p) for z' 
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configurations witti density p relax to tlie equilibrium con- 
figuration. To elaborate, we measure the relaxation time 
tq when H{cr) becomes as a function of p. As shown in 
figure m Tq increases when p is increased. From the graph, 
we can claim that Pc > 0.69. Due to the rapid increase in 
To at high densities, a better numerical estimation of the 
lower bound becomes difficult. 

Smaller upper bounds. The transition point is ob- 
tained from the thermodynamic limit of the probability of 
finding frozen particles in initial configurations. However, 
as discussed in the previous subsection, direct numerical 
estimation is practically impossible because the relaxation 
time becomes too long to distinguish between a slow relax- 
ation process and a freezing state. Nevertheless, in order 
to know the upper bounds of the transition point, we sim- 
ply require to observe a percolation of frozen particles in 
the direction from left to right under the boundary con- 
dition that particles are filled in the left- and right-side 
regions, while no particles exist in the top and bottom re- 
gions [22j. In fact, if such a directed percolation occurs in 
this modified system, frozen particles would be observed 
in the original system with the filled boundary conditions. 

In order to detect a directed percolation, we numerically 
investigate the probability Pl{p) of finding a frozen con- 
figuration in the modified system with M — L^/^, where 
< z""'^ < 1 because M < L. Then, for an assumed value 
of z, we determine whether a density pdp exists such that 
limi_s.oo Pi(yo) = 9{p — /Cdp), where 9{ ) is the Heaviside 
function. If such a density pdp exists, it corresponds to a 
directed percolation point with the assumed value z. In 
this manner, we can determine a set of values (pdp,2), 
and pdp provides an upper bound of the transition point. 
Practically, in numerical experiments, for a given value 
of z, we consider a cross point of the two curves Pl/2{p) 
and Pl{p), which is denoted by p^j.{z~^). We then check 
whether or not p^j.{z^^) pdp in the limit i — > oo. 

We first set z~-^ = 0.63, which is the dynamical expo- 
nent of the standard directed percolation [21]. As shown 
in figure [31 the numerical data up to i = 4000 suggest 
that p^i- are too scattered to judge convergence to a spe- 
cial value where the density is changed by an increment of 
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Fig. 4; Pl(p) with z^'^ = 0.64 and 0.65 (top), z~' 
0.67 (middle), and z^'^ = 0.68 and 0.69 (bottom). 



0.66 and 



0.0005. We thus seek the possibility of the convergence of 

^. Specifically, by changing 



p^j. for a different value of z 



z 



-1 



with an increment of 0.01 from 0.63, we measured 



J2i=iiP"' ~ Pci)^/3 for each value 
Here, we estimated p^/ by approximating Pl{p) 



Pl{p)- The resuhs for 0.64 < < 0.69 are shown 
in figure [31 From the graphs, one may conjecture, for 
example, that a good convergence of p^j. is expected for 
z~^ = 0.68. We then quantitatively investigate the extent 
of the convergence b^ measuring the average value of the 
cross points pcr = Per with Li — 500 x 2* and 

the deviation Apc 
of 2 

as a piece-wise linear function. When Apcr is sufficiently 
small, we expect an obvious convergence of p^j.{z~^) in the 
limit L ^ oo. 

Figure [3 presents the numerical results of Apcr for 
0.63 < < 0.78, where Apcr exhibits an oscillatory be- 
haviour. Such a behaviour has never been reported. Here, 
we assume that the local minima in Apcr provide con- 
vergence points, by relying on the interpretation that the 
existence of such multi-convergence points implies the co- 
existence of different types of directed percolations with 
different exponents z at different densities pdp- This 
leads to the conclusion that the convergence points are 
z~^ = 0.64,0.68,0.73, and 0.75, within numerical accu- 
racy. The corresponding values of pdp are displayed in 
figure [6l We also conjecture that there are more conver- 
gence points at a larger value of z~^, even though it is not 
easy to obtain them numerically due to the longer numer- 
ical simulations required. In principle, there exists a least 
upper bound pdp*, that can be obtained by this procedure. 
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Fig. 5: Apcr for z ^ . 

and a corresponding dynamical exponent . From numer- 
ical experiments, we have pdp* < 0.72f7 and z* < 1/0.75. 
This result strongly suggests that the freezing transition 
in our model is not characterized by a combination of ex- 
ponents of a standard directed percolation, such as the 
dynamical exponent Zdp — 1/0.63, unlike in the case of 
the TBF class il^IIHHinilll] ■ 

Concluding remarks. — The main achievement of 
this work is the presentation of a direction-free KCM that 
exhibits a freezing transition on the square lattice. We 
conjecture by numerical experiments that the singular be- 
haviour at the freezing transition does not belong to the 
TBF class P^ITMSnil^ . Lastly, we comment on future 
studies in the following text. 

First, since we have focused on the limiting case T — 0, 
it is not clear what happens in finite temperature cases. 
With regard to this problem, we note that the transition 
point above which frozen particles appear is independent 
of temperature because it is determined for each initial 
configuration, while the equilibrium density depends on 
temperature. Then, if frozen particles are not involved in 
an equilibrium configuration, the system might not reach 
the equilibrium configuration and remain trapped in a 
freezing state. Here, numerical simulations indicate that 
the fraction of frozen particles changes discontinuously at 
the transition point and that the fraction above the tran- 
sition point is always more than 0.5, which is an upper 
bound of equilibrium densities for any T . Therefore, we 
conjecture that a freezing transition occurs at the transi- 
tion point independent of the temperature. This conjec- 
ture will be confirmed in a future study. 

Second, the characterization of the singular behaviour 
near the freezing transition in our model remains to be 
solved. The correct value or more precise estimations for 
the transition point pc and the dynamical exponent z need 
to be derived theoretically and numerically. In particular, 
it might be interesting if one explicitly constructs a di- 
rected percolation problem related to the behavior near 
the transition point. Furthermore, the manner of diver- 
gence of the relaxation time tq should be clarified. We 
conjecture that tq in our model exhibits a Vogel-Fulcher 



Fig. 6: Dependence of estimated pdp for z ^ = 0.64, 0.68, 0.73, 
and 0.75. 

type singularity when p approaches from below, in a 
manner similar to that in the TBF class. Presently, we do 
not have clear evidence for the conjecture, because pc has 
not been estimated with sufficient accuracy as yet. After 
obtaining precise estimations for p^ and z, we will be able 
to validate our conjecture by numerical experiments. 

Third, the theoretical analysis of the model on a Bethe 
lattice will be performed in order to enhance the un- 
derstanding of the universality of freezing transitions in 
KCMs. A concrete question is whether the singular be- 
haviour of the freezing transitions observed here can be 
characterized by power-law exponents associated with a 
mode-coupling equation, as discussed in the Fredrickson- 
Andersen model [T5ll23l - [27] . If the answer is yes, since it is 
different from the Vogel-Fulcher type singularity, we will 
clarify the origin of the difference between the behaviours 
of the model on the Bethe lattice and on the square lattice. 

By addressing the points above, we wish to understand 
how freezing transitions in KCMs are related, or unre- 
lated, to jamming transitions and also wish to establish 
the universality classes associated with freezing transitions 
in non-equilibrium phenomena. 

if. if. 
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